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I.  INTRODUCTION 

In  order  to  investigate  the  radiation  patterns  of  antennas  in 
a  complex  environment  such  as  on  an  aircraft  the  following  Fortran 
IV  computer  code  has  been  developed.  The  computer  code  is  used  to 
compute  the  near  zone  radiated  fields  for  antennas  mounted  on  an  el¬ 
liptic  cylinder  and  in  the  presence  of  a  set  of  finite  flat  plates. 

The  analysis  applied  in  the  development  of  the  code  is  based  on  the 
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geometrical  theory  of  diffraction  (GTD)  ’  . 

The  code  allows  the  user  to  simulate  a  wide  variety  of  complex 
electromagnetic  radiation  problems  using  the  cylinder/plates  modal. 

For  example,  the  elliptic  cylinder  can  be  used  to  simulate  the  fuselage 
or  jet  intake  of  an  aircraft;  whereas,  the  plates  are  used  to  represent 
the  wings,  stabilizers,  stores,  etc.  Alternatively,  the  antenna  could 
be  mounted  directly  on  a  ship  mast.  In  this  case  the  mast  could  be 
approximated  by  the  elliptic  cylinder  with  the  other  ship  structures 
simulated  by  flat  plates.  Note  that  the  plates  can  be  attached  to 
the  cylinder  and/or  to  other  plates.  In  fact,  the  plates  can  be  connected 
together  to  form  a  box.  In  terms  of  special  sections  in  the  input 
data  set,  antennas  mounted  on  missile  configurations  will  be  discussed. 

The  T-tail  effect  and  mutual  coupling  effect  of  antenna  arrays  mounted 
on  nearly  flat  fuselage  are  also  illustrated  by  examples. 

As  with  any  ray  optical  solution  such  as  this  GTD  code,  there 
is  a  limit  to  the  number  of  interactions  included  in  the  field  computation. 
In  this  case,  the  code  includes  the  source,  reflected,  diffracted, 
reflected/reflected,  reflected/diffracted,  and  diffracted/reflected. 

In  addition,  the  user  can  request  the  code  to  add  various  diffracted/ 
diffracted  terms.  This  implies  that  the  code  can  handle  structures 
for  which  the  energy  does  not  significantly  bounce  back-and-forth  across 
the  target.  In  any  event,  the  code  automatically  shadows  all  terms, 
such  that  if  a  higher-order  interaction  should  have  been  included  the 
resulting  pattern  will  contain  a  discontinuity.  These  higher-order 
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terms  are  normally  negligible  and  can  only  affect  the  pattern  in  rather 
small  sectors.  However,  if  they  are  significant  in  some  region,  the 
amplitude  of  the  jump  is  associated  with  the  radiation  level  of  the 
missing  higher-order  term.  Consequently  when  the  solution  fails  because 
of  a  lack  of  higher-order  terms,  it  tends  to  Indicate  its  failure. 

The  code  has  the  flexibility  to  handle  arbitrary  pattern  cuts. 

In  addition,  an  arbitrary  antenna  type  can  be  analyzed  provided  the 
current  distribution  across  the  aperture  is  known.  This  is  done  by 
approximating  the  distribution  by  a  set  of  magnetic  current  elements 
mounted  on  the  cylinder  surface.  The  magnetic  current  elements  have 
a  cosine  distribution  along  the  magnetic  current  direction  and  are 
uniform  in  the  orthogonal  direction.  The  code  can,  also,  treat  a  mono¬ 
pole  or  monopole  arrays;  however,  length  of  each  element  cannot  exceed 
a  quarter  wavelength. 

The  mutual  coupling  effect  for  monopole  arrays  mounted  on  fuselage 
can  be  handled  by  thin-wire  theory^,  if  the  region  near  the  array  is 
nearly  flat.  For  engineering  purposes,  image  theory 
can  be  applied  to  calculate  the  relative  current  distributions  as 
equivalent  dipole  arrays.  The  relative  current  value  on  each  dipole 
Is  then  taken  to  be  part  of  the  input  data  for  each  monopole  source 
specification.  The  final  pattern  is  the  superposition  of  the  contri¬ 
butions  from  each  individual  monopole. 

The  limitations  associated  with  the  computer  code  result  from 
the  basic  nature  of  the  analyses.  The  solution  is  derived  using  the 
6TD,  which  is  a  high  frequency  approach.  In  terms  of  the  scattering 
from  plate  structures  this  means  that  each  plate  should  have  edges 
at  least  a  wavelength  long.  In  terms  of  the  cylinder  structure  its 
major  and  minor  radii  should  be  at  least  a  wavelength  in  extent.  In 
addition,  each  antenna  element  should  be  at  least  a  wavelength  from 
all  edges.  In  many  cases,  the  wavelength  limit  can  be  reduced  to  a 
quarter  wavelength  for  engineering  purposes. 
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The  present  form  of  the  computer  code  is  not  large  in  terms  of 
computer  storage  and  executes  a  pattern  for  a  single  antenna  element 
very  efficiently.  The  present  code  requires  approximately  200  K  bytes 
of  storage.  It  will  run  a  pattern  cut  of  360  points  for  a  simple  air¬ 
craft  model  with  one  antenna  element  in  approximately  30  seconds  on 
a  CDC-6600  computer. 

This  user's  manual  is  designed  to  give  an  overall  view  of  the 
operation  of  the  computer  code,  to  instruct  a  user  in  how  to  use  it 
to  model  structures,  and  to  show  the  validity  of  the  code  by  comparing 
various  computed  results  against  measured  data  whenever  available. 

Section  II  describes  an  overall  view  of  the  organization  of  the  program. 
The  definition  of  the  input  is  given  in  Section  III.  How  to  apply 
the  capabilities  of  this  input  data  to  a  practical  structure  is  briefly 
discussed  in  Section  IV.  This  includes  a  clarification  of  the  subtle 
points  of  interpreting  the  input  data.  The  representation  of  the  output 
is  discussed  in  Section  V.  Various  sample  problems  are  presented  in 
Section  VI  and  VII  to  illustrate  the  operation,  versatility,  and  validity 
of  the  code. 


II.  PRINCIPLES  OF  OPERATION 

The  analytical  modeling  of  complex  scattering  shapes  in  order 

to  predict  the  radiation  patterns  of  antennas  has  been  accomplished 
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by  the  use  of  the  Geometrical  Theory  of  Diffraction  (GTD)  *  .  This 
is  a  high  frequency  technique  that  allows  a  complicated  structure  to 
be  approximated  by  basic  shapes  representing  canonical  problems  in 
the  GTD.  These  shapes  include  flat  and  curved  wedges  and  convex  curved 
surfaces.  The  GTD  is  a  ray  optical  technique  and  it  therefore  allows 
one  to  gain  some  physical  insight  into  the  various  scattering  and  dif¬ 
fraction  mechanisms  involved.  Consequently,  one  is  able  to  quickly 
seek  out  the  dominant  or  significant  scattering  and  diffraction  mechanism 
for  a  given  geometrical  configuration.  This,  in  turn,  leads  to  an 
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accurate  engineering  solution  to  practical  antenna  problems.  This 
approach  has  been  used  successfully  in  the  past  to  model  aircraft 
shapes'*’4’*’’®  and  ship-like  structures7’®*®. 

This  section  briefly  describes  the  basic  operation  of  this  code 
for  the  analysis  of  airborne  antenna  patterns.  The  present  version 
of  the  code  allows  analysis  of  structures  that  can  be  modeled  by  con¬ 
cave  flat  plate  structures,  and  an  elliptic  cylinder  all  of  which  are 
built  up  from  the  basic  canonical  structures.  These  shapes  allow  one 
to  model  a  wide  variety  of  structures  in  the  UHF  range  and  above  where 
the  scattering  structures  are  large  in  terms  of  a  wavelength.  The 
general  rule  is  that  the  lower  frequency  limit  of  this  solution  is 
dictated  by  the  spacings  between  the  various  scattering  centers  and 
their  overall  size.  In  practice  this  means  that  the  smallest  dimen¬ 
sions  should  be  on  the  order  of  a  wavelength.  This  can  often  be  relaxed 
to  approximately  a  quarter- wave length. 

The  positive  time  convention  eJa)t  has  been  used  in  this  scattering 
code.  Also,  the  radiation  patterns  are  assumed  to  be  in  the  near  field. 
All  of  the  structures  in  this  code  are  assumed  to  be  perfectly  conducting 
and  the  surrounding  medium  is  free  space. 

As  mentioned  above,  the  GTD  approach  is  ideal  for  a  general  high 
frequency  study  of  aircraft  antennas  in  that  only  the  most  basic 
structural  features  of  an  otherwise  very  complicated  structure  need 
to  be  modeled.  This  is  because  ray  optical  techniques  are  used  to 
determine  components  of  the  field  incident  on  and  diffracted  by  various 
structures.  Components  of  the  diffracted  fields  are  found  using  the 
GTD  solutions  in  terms  of  the  individual  rays  which  are  summed  with 
the  geometrical  optics  terms  at  the  field  point.  The  rays  from  a  given 
scatterer  tend  to  interact  with  other  structures  causing  various  higher- 
order  terms.  In  this  way  one  can  trace  out  the  various  possible  com¬ 
binations  of  rays  that  interact  between  scatterers  and  determine  and 
Include  only  the  dominant  terms.  Thus,  one  need  only  be  concerned 
with  the  important  scattering  components  and  neglect  all  other  higher- 
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order  terms.  This  method  leads  to  accurate  and  efficient  computer 
codes  that  can  be  systematically  written  and  tested. 

Complex  problems  are  built  up  from  similar  components  in  terms 
of  a  modular  computer  code.  This  modular  approach  is  illustrated 
in  the  block  diagram  of  the  main  program  shown  in  Table  I.  The  code 
is  broken  up  into  many  subroutines  that  represent  different  scattered 
field  components,  ray  tracing  methods,  and  shadowing  routines.  As 
can  be  seen  from  the  flow  chart,  the  code  is  structured  so  that  all 
of  one  type  of  scattered  field  is  computed  at  one  time  for  the  complete 
pattern  cut  so  that  the  amount  of  core  swapping  is  minimized  thereby 
reducing  overlaying  and  increasing  efficiency.  This  also  is  an  important 
feature  that  allows  the  code  to  be  used  on  small  computers  that  are 
not  large  enough  to  accept  the  entire  code  at  one  time.  The  code  can 
be  broken  into  smaller  overlay  segments  which  will  individually  fit 
in  the  machine.  The  results  are,  then,  superimposed  in  the  main  program 
as  the  various  segments  are  executed. 

The  subroutines  for  each  of  the  scattered  field  components  are 
all  structured  in  the  same  basic  way.  First,  the  ray  path  is  traced 
backward  from  the  chosen  observation  point  to  a  particular  scatterer 
and  subsequently  to  the  source  using  either  the  laws  of  reflection 
or  diffraction.  Each  ray  path,  assuming  one  is  possible,  is  then  checked 
to  see  if  it  is  shadowed  by  any  structure  along  the  complete  ray  path. 

If  it  is  shadowed  the  field  is  not  computed  and  the  code  proceeds  to 
the  next  scatterer  or  observation  point.  If  the  path  is  not  interrupted 
the  scattered  field  is  computed  using  the  appropriate  GTD  solutions. 

The  fields  are  then  superimposed  in  the  main  program.  This  shadowing 
process  is  often  speeded  up  by  making  various  decisions  based  on  bounds 
associated  with  the  geometry  of  the  structure.  This  type  of  knowledge 
is  used  wherever  possible. 
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The  shadowing  of  rays  is  a  very  important  part  of  the  GTD  Scat¬ 
tering  code.  It  is  obvious  that  this  approach  should  lead  to  various 
discontinuities  in  the  resulting  pattern.  However,  the  GTD  diffraction 
coefficients  are  designed  to  smooth  out  the  discontinuities  in  the 
fields  such  that  a  continuous  field  is  obtained.  When  a  scattered 
field  is  not  included  in  the  result,  the  lack  of  its  presence  is  apparent. 
This  can  be  used  to  advantage  in  analyzing  complicated  problems. 

Obviously  in  a  complex  problem  not  all  the  possible  scattered  fields 
can  be  included.  In  the  GTD  code  the  importance  of  the  neglected  terms 
are  determined  by  the  size  of  the  so-called  gliches  or  jumps  in  the 
pattern  trace.  If  the  gliches  are  small  no  additional  terms  are  needed 
for  a  good  engineering  solution.  If  the  gliches  are  large  it  may  be 
necessary  to  include  more  terms  in  the  solution.  In  any  case  the  user 
has  a  gauge  with  which  he  can  examine  the  accuracy  of  the  results  and 
is  not  falsely  led  into  believing  a  result  is  correct  when  in  fact 
there  could  be  an  error. 

The  brief  discussion  of  the  operation  of  this  code  given  above 
should  help  the  user  get  a  feel  for  the  overall  code  so  he  might  better 
understand  the  code's  capabilities  and  interpret  its  results.  The 
code  is  designed,  however,  so  that  the  general  user  can  run  the  code 
without  knowing  all  the  details  of  its  operation.  Yet,  he  must  become 
familiar  with  the  input/output  details  which  will  be  discussed  in  the 
next  three  sections. 
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TABLE  I 

BLOCK  DIAGRAM  OF  THE  MAIN  PROGRAM 

(start) 

SET  DEFAULT  DATA 

I 

READ  INPUT  DATA 
— ►  ECHO  INPUT  DATA 

(SEE  TABLE  II) 

1  IF  READ  INPUT  COMPLETE  1  YES 


NO 


CONVERT  UNITS  TO  WAVELENGTHS 
DEFINE  PATTERN  COORDINATE  INFORMATION 
INITIALIZE  TOTAL  FIELD  TO  ZERO 


DO  LOOP  ON  ANTENNA  ELEMENTS 


SPECIFY  ANTENNA  GEOMETRY  AND  TYPE 
DETERMINE  ALL  FIXED  GEOMETRY 


COMPUTE 

VARIOUS 

GTD  TERMS  1 

(NOTE: 

PATTERN  LOOP  IN  EACH  GTD  TERM  SUBROUTINE) 

a) 

SOURCE 

b) 

REFLECTED 

c) 

DIFFRACTED 

d) 

REFLECTED/REFLECTED 

e) 

REFLECTED/DIFFRACTED 

f) 

DIFFRACTED/REFLECTED 

COMPUTED  SPECIFIED  DOUBLE  DIFFRACTION 


COMPUTED  SPECIFIED  LINE  SOURCE  ARRAY 


CONVERT  X, Y,Z  FIELD  COMPONENTS  TO 
THETA  AND  PHI  IN  PATTERN 
COORDINATE  SYSTEM 


PRINT  AND/OR  PLOT  THETA  AND  PHI 
FIELD  COMPONENTS 
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III.  DEFINITION  OF  INPUT  DATA 

The  method  used  to  input  data  into  the  computer  code  is  presently 
based  on  a  command  word  system.  This  is  especially  convenient  when 
more  than  one  problem  is  to  be  analyzed  during  a  computer  run.  The 
code  stores  the  previous  input  data  such  that  one  need  only  input  that 
data  which  needs  to  be  changed  from  the  previous  execution.  Also, 
there  is  a  default  list  of  data  so  for  any  given  problem  the  amount 
of  data  that  needs  to  be  input  has  been  shortened.  The  organization 
of  the  input  data  is  illustrated  in  Table  II. 

In  this  system,  all  linear  dimensions  may  be  specified  in  either 
meters,  inches,  or  feet  and  all  angular  dimensions  are  in  degrees. 

All  the  dimensions  are  eventually  referred  to  a  fixed  cartesian  coordinate 
system  used  as  a  common  reference  for  the  source  and  structures.  There 
is,  however,  a  geometry  definition  coordinate  system  that  may  be  defined 
using  the  "RT:"  command.  This  command  enables  the  user  to  rotate  and 
translate  the  coordinate  system  to  be  used  to  input  any  selected  data 
set  into  the  best  coordinate  system  for  that  particular  geometry. 

Once  the  "RT:"  command  is  used  all  the  input  following  the  command 
will  be  in  that  rotated  and  translated  coordinate  system  until  the 
"RT:"  command  is  called  again.  The  only  exception  to  this  is  that 
the  cylinder  will  always  be  in  the  reference  coordinate  system.  See 
below  for  more  details.  There  is  also  a  separate  coordinate  system 
that  can  be  used  to  define  a  pattern  coordinate  system.  This  is  discussed 
in  more  detail  in  Section  III-C  in  terms  of  the  "PD:"  command. 

It  is  felt  that  the  maximum  usefulness  of  the  computer  code  can 
be  achieved  using  it  on  an  interactive  computer  system.  As  a  consequence, 
all  input  data  are  defined  in  free  format  such  that  the  operator  need 
only  put  commas  between  the  various  inputs.  This  allows  the  user  on 
'an  interactive  terminal  to  avoid  the  problems  associated  with  typing 
in  the  field  length  associated  with  a  fixed  format.  This  method  also 
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is  useful  on  batch  processing  computers.  Note  that  all  read  statements 
are  made  on  unit  #5,  i.e.,  READ  (5,*),  where  the  symbol  refers 
to  free  format.  Other  machines,  however,  may  have  different  symbols 
representing  free  format. 

In  all  the  following  discussions  associated  with  logical  variables 
a  "T"  will  imply  true,  and  an  "FM  will  imply  false.  The  complete  words 
true  and  false  need  not  be  input  since  most  compilers  just  consider 
the  first  character  in  determining  the  state  of  the  logical  variable. 

The  following  list  defines  in  detail  each  command  word  and  the 
variables  associated  with  them.  Section  VI  will  give  specific  examples 
using  this  input  method.  Note  that  the  program  halts  execution  by 
sensing  the  end-of-file  mark  associated  with  the  input  data  set. 


TABLE  II 


BLOCK  DIAGRAM  OF  THE  INPUT  DATA 
ORGANIZATION  FOR  THE  COMPUTER  CODE 


(start) 


w 


Initialize  Default  Data 

I 


Read  and  Write  Command  Word 


7 


MM 


NADC-A0096-30 
A.  Command  TO: 


This  command  enables  the  user  to  obtain  an  extended  output  of 
various  intermediate  quantities  in  the  computer  code.  This  is  useful 
in  testing  the  program  or  in  analyzing  the  contributions  from  various 
scattering  mechanisms  in  terms  of  the  total  solution. 

1.  READ:  LDEBUG,  LTEST,  LOUT 

a)  LDEBUG:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  debug  the  program  if  errors  are 
suspected  within  the  program.  If  set  true,  the 
program  prints  out  data  on  unit  #6  associated 
with  each  of  its  internal  operations.  These 
data  can,  then,  be  compared  with  previous  data 
which  are  known  to  be  correct.  It  is,  also. 
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used  to  insure  initial  operation  of  the  code. 

Only  one  pattern  angle  is  considered. 

(normally  set  false) 

b)  LTEST:  This  a  logical  variable  defined  by  T  or  F.  It 

is  used  to  test  the  input/output  associated 
with  each  subroutine.  The  data  written  out 
on  unit  #6  are  associated  with  the  data  in  the 
window  of  the  subroutine.  They  are  written 
out  each  time  the  subroutine  is  called.  It 
is,  also,  used  to  insure  initial  operation  of 
the  code.  Only  one  pattern  angle  is  considered, 
(normally  set  false) 

c)  LOUT:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  output  data  on  unit  #6  associated 
with  the  main  program.  It, also,  is  used  to  ini¬ 
tially  insure  proper  operation.  It  can  be  used 
to  examine  the  various  components  of  the  pattern, 
(normally  set  false) 

READ:  LTERM,  LCORNR 

a)  LTERM:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  tell  the  code  whether  or  not  in¬ 
dividual  terms  are  desired  during  the  compu¬ 
tation. 

(normally  set  false) 

b)  LCORNR:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  tell  the  code  whether  or  not  corner 
diffraction  is  desired  during  the  computation, 
(normally  set  true) 


14 


NADC-80096-30 


3.  READ:  (LTRM(J),  J=l,7) 

a)  LTRM{J):These  are  logical  variables  defined  by  T  or 
F  to  specify  a  set  of  individual  scattering 
components  that  are  to  be  inlcuded  in  the  scat¬ 
tered  field  computation.  The  components  are 
defined  by  the  following  number  designations. 

J=l:  source  field 
J=2:  single  reflected  field 
J=3:  single  diffracted  field 
J=4:  double  reflected  field 
J=5:  ref lected-diffracted  field 
J=6:  diffracted-ref lected  field 
J=7:  double  diffracted  field 

4.  READ:  MPI,  MPF,  MPS 

a)  MPI, MPF, MPS:  These  are  integer  variables  to  define 
plates  used  in  computation,  where 
MPI  =  initial  plate 
MPF  =  final  plate 

MPS  =  increment  in  plates  going  from  initial 
to  final  plate 

Note:  MPI=1,  MPF=3,  and  MPS=2  implies  plates  1  and  3  are 
included  in  the  computation. 
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5.  READ:  (MEI(MP),  MEF(MP),  MES(MP),  MP1 * 3MPI,MPF,MPS) 

a)  MEI(MP),MEF(MP),MES(MP):  These  are  dimensioned  integer 
variables  to  define  the  edges  on  the  MPth  plate 
used  in  the  computation,  where 
MEI(MP)  =  Initial  edge  on  plate  MP 
MEF(MP)  3  final  edge  on  plate  MP 
MES(MP)  =  increment  in  edges  going  from  MEI(MP) 
to  MEF(MP) 

B.  COMMAND  UN: 


This  command  enables  the  user  to  specify  the  units  used  for  all 
following  linear  dimensions  in  input  data  list. 

1.  READ:  I UNIT 

A)  IUNIT:  This  is  an  integer  variable  that  defines  the 
units.  If 

1  -►  meters 

IUNIT  =  ■  2  -  feet 

3  -►  inches 
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C.  COMMAND  PD: 


This  command  enables  the  user  to  define  the  pattern  axis  of  ro¬ 
tation,  the  angular  range,  range  from  origin  to  receiver,  and  the  fre¬ 
quency  for  the  desired  conical  pattern. 

This  set  of  data  is  associated  with  the  conical  pattern  desired 
during  execution  of  the  program.  The  pattern  axis  is  defined  by  the 
spherical  angles  (THC,  PHC)  as  illustrated  in  Figure  1.  These  angles 
define  a  radial  vector  direction  which  points  in  the  direction  of  the 
pattern  axis  of  rotation.  These  angles  actually  set-up  a  new  coordi¬ 
nate  system  in  relation  to  the  original  fixed  coordinates.  The  new 
cartesian  coordinates  defined  by  the  subscript  "p"  are  found  by  first 
rotating  about  the  z-axis  the  angle  PHC  and,  then,  about  the  yp-axis 
the  angle  THC.  The  pattern  is,  then,  taken  in  the  "p"  coordinate  sys¬ 
tem  in  terms  of  spherical  angles.  The  theta  angle  of  the  pattern  taken 
about  the  zp-axis  is  defined  by  THETA.  The  phi  angle  is  defined  by  the 
next  read  statement.  In  the  present  form  the  program  will,  then,  com¬ 
pute  any  conical  pattern  in  that  THETA  is  used  as  the  conical  pattern 
angle  about  the  zp-axis  for  the  complete  pattern  calculation. 

As  an  aid  in  setting  up  the  "p"  coordinate  system  the  following 
set  of  equations  give  the  relationships  between  (THC,  PHC)  and  the 
xp,  yp,  zp-axes.  Note  that  the  "p"  axes  are  defined  as  radial  vector 
directions  in  a  spherical  coordinate  system: 
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5  =  cos ( PHC ) s i n ( THC+900 ) x+ s i n ( PHC ) s i n ( THC+900 ) y+cos ( THC+900 ) z 

yp  =  cos(PHC+90°)x+sin(PHC+90°)y 

zp  =  cos(PHC)sin(THC)x+sin(PHC)sin(THC)y+cos(THC)z 

where  0  £  THC  £  180°  and  0  £  PHC  £  360°.  In  its  present  form  it  should 
be  noted  that  the  user  may  not  be  able  to  define  the  xp-axis  at  the 
starting  location  that  he  desired.  In  addition,  the  rotation  of  the 
pattern  may  be  in  the  opposite  sense  using  this  approach.  However, 
these  problems  can  be  easily  overcome  with  properly  written  plot  routines 


1.  READ:  THC,  PHC,  THETA 

a)  THC, PHC:  These  are  real  variables.  They  are  input 

in  degrees  and  define  the  axis  of  rotation 
about  which  a  conical  pattern  will  be  com¬ 
puted. 


b)  THETA:  This  is  a  real  variable.  It  is  input  in 

degrees  and  used  to  define  the  conical  angle 
about  the  axis  of  rotation  for  the  desired 
pattern. 


■m's 
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X 

Figure  1.  Definition  of  pattern  axis. 


2.  READ:  IPS,  IPF,  IPD 

a)  IPS, IPF, IPD:  These  are  integer  variables  used  to  define 
angles  in  degrees.  They  are,  respectively, 
the  beginning,  ending,  and  incremental 
values  of  the  phi  pattern  angle. 

As  a  result  of  the  input  given  by  the  two  previous  read 
statements,  the  operator  has  completely  defined  the  desired 
conical  pattern  to  be  computed  during  execution  of  the  pro¬ 
gram. 
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3.  READ:  R,  FRQG 

a)  R:  This  is  a  real  variable  which  is  used  to 

define  the  range  in  linear  units  from  the 
origin  to  the  reciever.  If  the  pattern  is 
to  be  taken  in  the  far  field  of  the  structure 
one  must  be  very  careful  how  he  inputs  the 
range.  This  point  is  discussed  in  detail 
in  the  next  section. 

b)  FRQG:  This  is  a  real  variable  which  is  used  to 

define  the  frequency  in  gigahertz. 

D.  COMMAND  RT: 


This  command  enables  the  user  to  translate  and/or  rotate  the 
coordinate  system  used  to  define  the  input  data  in  order  to  simplify 
the  specif ication  of  the  plate  geometry.  The  geometry  is  illustrated 
in  Figure  2. 
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Figure  2.  Definition  of  rotate-translate  coordinate  system  geometry. 

1.  READ:  (TR(N) ,  N=l,3) 

a)  TR(N):  This  is  a  dimensioned  real  variable.  It 

is  used  to  specify  the  origin  of  the  new 
coordinate  system  to  be  used  to  input  the 
data  for  the  plate  structures.  It  is  input 
on  a  single  line  with  the  real  numbers  being 
the  x,y, z  coordinates  of  the  new  origin  which 
corresponds  to  N=l,2,3  respectively. 


NADC -80096-30 


2.  READ:  THZR,  PHZR,  THXR,  PHXR 

a)  THZR, PHZR:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  zR-axis  of  the  new  coordinate  system  as 
if  it  was  a  radial  vector  in  the  reference 
coordinate  system. 

b)  THXR, PHXR:  These  are  real  variables.  They  are  input 

in  degrees  as  spherical  angles  that  define 
the  xR-axis  of  the  new  coordinate  system  as 
if  it  was  a  radial  vector  in  the  reference 
coordinate  system. 


The  new  xR-axis  and  zR-axis  must  be  defined  orthogonal  to  each 
other.  The  new  yR-axis  is  found  from  the  cross  product  of  the  Xp  and 
Zp  axes.  All  the  subsequent  inputs  will  be  made  relative  to  this  new 
coordinate  system,  which  is  shown  as  (xR,  yR,  zR)  unless  command  "RT:" 
is  called  again  and  redefined. 


E. 


COMMAND  PG: 


m 
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This  command  enables  the  user  to  define  the  geometry  of  the  flat 
plate  structures  to  be  considered.  The  geometry  is  illustrated  in 
Figure  3.  It  can  be  called  repeatedly  up  to  25  times. 


Figure  3.  Definition  of  flat  plate  geometry. 

1.  READ:  MCX(MPX),  LATACH(MPX) 

a)  MCX(MPX):  This  is  a  dimensioned  integer  variable.  It 
is  used  to  define  the  number  of  corners  (or 
edges)  on  the  MPXth  plate.  Presently, 
1<MCX(MPX)<6  with  1<MPX<25 
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b)  LATACH(MPX):  This  is  a  logical  variable  defined  by 
T  or  F.  It  is  used  to  indicate  if  the 
MPXth  plate  is  attached  to  the  fuselage 
(T)  or  not  (F).  Note  that  all  attached 
plates  should  be  defined  within  the 
first  six  plates.  The  first  and  last 
corners  of  attached  plates  should  be 
specified  on  or  near  the  fuselage. 

2.  READ:  (PVC(N,ME,MPX),N=1,3) 

As  stated  earlier  the  locations  of  the  corners  of  the  flat 
plates  are  input  in  terms  of  the  x,  y,  z  coordinates  in 
the  specified  cartesian  coordinate  system. 

a)  PVC(N,ME,MPX):  This  is  a  triply  dimensioned  real  vari¬ 
able.  It  is  used  to  specify  the  location 
of  the  MEth  corner  of  the  MPXth  plate. 

It  is  input  on  a  single  line  with  the 
real  numbers  being  the  x,y,z  coordinates 
of  the  corner  which  correspond  to  N=l, 
2,3,  respectively,  in  the  array.  For 
example,  the  array  will  contain  the 
following  for  plate  #1  and  corner  #2 
located  at  x=2.,  y=4. ,  z=6.: 

PVC(1,2,1)  *  2. 

PVC (2,2,1)  =  4. 

PVC(3,2,1)  =  6. 

This  data  is  input  as:  2. ,4. ,6. 


24 


NADC-80096-30 


Considering  the  flat  plate  structure  given  in  Figure  3,  the  input 
data  is  given  by 


1.,  1.,  0 

:  corner  #1 

-1.,  1.,  0. 

:  corner  #2 

-l.,-l.,  0. 

:  corner  #3 

l.,-l.,  0. 

:  corner  #4 

Presently: 

1  <  MPX  <25 

1  <  ME  <  6 

i  <  N  <3 

plate  #1 


(See  Chapter  IV  for  further  details  in  defining  the  corner  points.) 


F. 


COMMAND  SG: 


1 


/  READ:  MSX  7 

_  € — 


DO  3402  MS*1,MSX 

—  i  — 7 

/READ:  PHSA(MS) ,  ZSA(MS)  / 

EZ 

READ:  SLOTAA(MS),  SLOTBA(MS) , 
BETADA(MS),  SMONOA(MS), 
(JANTA{NA,MS),NA=1,NAX) 

/  READ:  WMA(MS),  WPA(MS)  / 

- "rzz — 


J 3402  CONTINUE 


This  command  enables  the  user  to  specify  the  location  and  type 
of  antenna  to  be  used.  The  geometry  is  Illustrated  in  Figure  4. 
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Figure  4.  Definition  of  the  source  location  for  computer  code. 
Note  that  PHSA  =  $  and  ZSA  =  -  |Z  |  in  the  above 
drawings. 
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1.  READ:  MSX 

a)  MSX:  This  is  an  integer  variable  which  defines 

the  maximum  number  of  elemental  radiators 
to  be  considered  during  execution  of  the 
program.  Presently,  1  5  MSX  <  10. 

2.  READ:  PHSA(MS),  ZSA(MS) 

a)  PHSA(MS),ZSA(MS):  These  are  real  variables  used  to 
specify  the  phi -angle  (in  degrees)  and  2 
location  of  MSth  antenna.  (Refer  to  Fig¬ 
ure  4)  Note:  -90°  <  PHSA  <  90°. 

3.  READ:  SLOTAA(MS),  SLOTBA(MS),  BETADA(MS) ,  SMONOA(MS), 

(JANTA(NA.MS),  NA=1,NAX) 

a)  SLOTAA(MS),  SLOTBA(MS):  These  are  real  variables  used 

to  specify  the  narrow  (parallel  with  E  field) 
and  broad  (perpendicular  to  E  field)  dimensions 
of  the  slot  in  wavelengths. 

b)  BETADA(MS):  This  is  a  real  variable  used  to  specify 

the  angle  (in  degrees)  of  the  slot  relative 
to  the  fuselage  axis.  If  BETADA=0.,  then  is 
an  axial  slot.  If  BETADA=90.,  then  is  a 
circumferential  slot. 

c)  SMONOA(MS):  This  is  a  real  variable  used  to  specify 

the  length  of  monopole  in  wavelengths.  Note 
that  SMONOA  should  not  exceed  a  quarter  wave¬ 
length. 
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d)  JANTA{1,MS):  This  is  an  integer  variable  used  to  specify 

the  type  of  antenna  considered  in  computation: 

1  -►  arbitrarily  oriented  slot 
JANTA  =  <  2  -*■  circumferential  slot 
3  -*•  radial  monopole. 

4.  READ:  WMA(MS),  WPA(MS) 

a)  WMA(MS),  WPA(MS):  These  are  real  variables  used  to  specify 
the  magnitude  and  phase  (in  degrees)  of  excitation 
on  the  MSth  antenna.  If  an  array  is  used,  then 
the  excitation  including  the  coupling  effect 
on  the  radiators  may  be  obtained  using  a  thin- 
wire  code  as  shown  in  the  results  section. 

G.  COMMAND  RS: 

This  command  enables  the  user  to  reset  the  input  data  to  the 
default  case.  There  is  no  input  data  associated  with  this  command. 
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1.  READ:  N00X 

a)  NOOX:  This  Is  an  Integer  variable  that  specifies 

the  total  number  of  double  diffraction  terms 
desired.  Presently,  0^IDDX<10. 

2.  READ:  MPIDD(NOD),  MEIDD(NDD),  MP2DD(NDD),  ME2DD(NDD) 

a)  MPIOO(NDD),  MEIDD(NDD):  These  are  integer  dimensioned 
arrays  used  to  specify  the  plate  and  edge 
number  from  which  the  first  diffraction  oc¬ 
curs. 


b)  MP2DD(NDD),  ME2DD(NDD):  These  are  integer  dimensioned 
array  used  to  specify  the  plate  and  edge 
number  from  which  a  second  diffraction  occurs. 


I.  COMMAND  F6: 


/  READ: 
/  READ: 


AF,BF,CF,DF 

1 


(PVO(H),N«l,3) 

CL) 


7 

7 


This  command  enables  the  user  to  model  the  fuselage  by  a  composite 
elliptic  cylinder. 
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1.  READ:  AF,  BF,  CF,  DF 

a)  AF,  BF,  CF,  DF:  These  are  real  variables  that  specify 
the  semi-major  and  semi -minor  axes  of  the 
composite  elliptic  cylinder  used  to  model 
the  fuselage  as  shown  in  Figure  5.  Note 
how  this  figure  compared  with  Figure  4. 


Figure  5.  Definition  of  fuselage  geometry. 


2.  READ:  (PVO(N),  N=l,3) 

(a)  PVO(N);  This  is  a  real  dimensioned  array  that  defines 
the  location  of  the  origin  about  which  the 
pattern  is  taken,  i.e.,  PVO(N)  =  (x,y,z). 

Note  that  the  conical  pattern  will  be  taken 
about  this  origin. 
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J.  COMMAND  IG: 


This  command  enables  the  user  to  define  the  jet  intake  geometry. 
The  geometry  is  illustrated  in  Figure  6. 


1.  READ:  C,D,AI,BI,NINT,LFUS 


a)  C:  This  is  a  real  variable  which  defines  the 

radius  of  the  fuselage. 

b)  D:  This  is  a  real  variable  which  defines  the 

radial  distance  from  the  fuselage  origin 
to  the  intake  center. 


c)  AI,BI:  These  are  real  variables  that  define  the 

semi-major  and  semi -minor  dimensions  of  the 
jet  intakes. 


d)  NINT:  This  is  an  integer  variable  that  defines 

the  number  of  jet  intakes  used  in  computation. 
It  is  assumed  that  the  intakes  are  uniformly 
positioned  around  the  fuselage. 

e)  LFUS:  This  is  a  logical  variable  defined  by  T  or 

F  to  specify  whether  the  antenna  system  is 
mounted  on  the  fuselage  (T)  or  intake  (F). 
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This  command  enables  the  user  to  specify  a  line  source  distri¬ 
bution  along  the  z-axis.  It  is  used  in  an  array  pattern  multiplication 
analysis.  This  command  applies  only  when  one  has  a  uniform  geometry 
along  the  axis  of  the  fuselage. 


1.  READ:  LSTERM 


a)  LSTERM: 


LSTERM=0: 


LSTERM=1: 


LSTERM=2: 


This  is  an  integer  variable  that  indicates 
the  type  of  line  source  distribution  treated. 
The  current  distribution  and,  therefore, 
the  following  inputs  vary  according  to  the 
following  table. 


NINP  j<|>N 

I(z)  =  t  I  I|N  e  N  6(z-(N-S5)Az) 
N=1  " 

2.  READ:  NINP.DELZ 

3.  READ:  CURM.CPHAS 


-jkz  cos9f 


Id)  *  [e-otz+re-i6e-a<2L-z)]e'  s 

2.  READ:  SLENG , DBATT , GAMM , BETA , THSCAN 


I(z)  = 


!  .  ^  2(z-L/2) 


N  -jkz  cos0< 


>  e 


2:  READ":  SLENG,  NPOW.THSCAN 


LSTERM=3: 


LSTERM=4: 


I(z)  = 
2 


pAcl 


-jkz  cos0. 


cos  r;  j 

READ:  SLENG, CTERM, NPOW.THSCAN 


cos 


„z.  N  -JMp-p0) 

r)  e 


I(z)  ■  (■ 

where  p  -Vp2*!*-*,,)2 

2.  READ:  SLENG,  NPOW,  RHOP,  ZOP 
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The  input  data  is  interpreted  as  follows. 

a)  NINP:  This  is  an  integer  variable  that  defines 

the  number  of  current  samples. 

b)  DELZ:  This  is  a  real  variable  (Az)  that  defines 

the  current  sample  spacing  in  wavelengths. 

c)  CURM,  CPHAS:  These  are  real  variables  that  define  the 

magnitude  (|I|N)  and  phase  (<J>N)  of  the  cur¬ 
rent  elements. 

d)  SLENG:  This  is  a  real  variable  (L)  that  defines 

the  length  of  the  linear  array,  in  wavelengths. 

e)  DBATT:  This  is  a  real  variable  that  defines  the 

attenuation  (in  dB)  along  the  total  length 
(SLENG)  of  the  array.  Note  that  a  is  related 
to  DBATT. 

f)  GAMM,  BETA:  These  are  real  variables  (r  and  B)  that  define 

the  magnitude  and  phase  (in  degrees)  of  the 
reflection  coefficient  at  the  end  of  the 
traveling  wave  antenna  (LSTERM=1). 

g)  THSCAN:  This  is  a  real  variable  that  defines  the 

scan  angle  (in  degrees)  of  the  array. 

h)  NPOW:  This  is  an  integer  variable  (N)  that  defines 

the  exponent  in  the  previous  equations. 
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i)  CTERM:  This  is  a  real  variable  that  defines  the 

constant  (C)  in  the  previous  equations. 

j)  RHOP,  ZOP:  These  are  real  variables  that  define  the 

phase  distribution  across  an  aperture.  Note 
that  RHOP  and  ZOP  are  specified  in  wavelengths. 
In  terms  of  the  previous  definition  for  the 
case  (LSTERM=4)  RH0P=po  and  Z0P=ZQ. 

L.  COMMAND  SA: 


This  command  enables  the  user  to  get  a  quick  array  pattern  if 
the  structure  (geometry,  antenna  type,  and  location)  has  circumferential 
symmetry.  The  program  is  set  to  take  the  pattern  about  the  z-axis 
(longitudinal  axis)  in  one  degree  steps  in  this  case.  This  command 
must  be  used  in  conjunction  with  the  IG:  command.  Otherwise,  the 
program  aborts.  Presently,  this  command  applies  only  when  2  <  NINT 
<  4,  where  NINT  is  number  of  intakes  as  defined  in  the  "IG:"  command. 
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1.  READ:  NSYM 


a)  NSYM: 


This  is  an  integer  variable  which  defines 
the  number  of  arrays  placed  around  the  cir¬ 
cumference.  Presently,  2  <  NSYM  <  4.  It 
is  used  according  to  the  following  table: 


NSYM  =  | 


2  or  4  if  NINT  =  4 

3  if  NINT  =  3 

2  if  NINT  =  2 


2.  READ:  SAM, SAP 


a)  SAM,  SAP:  These  are  real  variables  used  to  specify 

the  magnitude  and  phase  (in  degrees)  of  the 
Nth  circumferential  array  element  excitation. 


M.  COMMAND  LP: 

This  command  enables  the  user  to  obtain  a  line  printer  listing 

of  the  total  fields  (Efl_,  E^).  No  input  is  required. 

yp  (pp 

N.  COMMAND  PP: 

9 

/  READ:  RADIUS,  IPLT  / 

^T_/ 

This  command  enables  the  user  to  obtain  a  pen  plot  of  the  total 

ffe'ds  <v  V- 
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1.  READ:  RADIUS,  IPLT 


a)  RADIUS:  This  is  a  real  variable  that  is  used  to  specify 

the  radius  of  the  polar  plot. 


b)  IPLT:  This  is  an  integer  variable  that  indicates 

the  type  of  polar  plot  desired,  such  that 


IPLT  = 


1  -*■  field  plot 

2  -►  power  plot 

3  >  dB  plot 


0.  COMMAND  EX: 


This  command  is  used  to  execute  the  code  so  that  the  total  fields 
may  be  computed.  After  execution  the  code  returns  for  another  possible 
command  word. 


This  concludes  the  definition  of  all  the  input  parameters  to 
the  program.  The  program  would,  then,  run  the  desired  data  and  output 
the  results  on  unit  #6.  However  as  with  any  sophisticated  program, 
the  definition  of  the  input  data  is  not  sufficient  for  one  to  fully 
understand  the  operation  of  the  code.  In  order  to  overcome  this  dif¬ 
ficulty  the  next  chapter  discusses  how  the  input  data  are  interpreted 
and  used  in  the  program. 
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IV.  INTERPRETATION  OF  INPUT  DATA 

This  computer  code  is  written  to  require  a  minimum  amount  of 
user  information  such  that  the  burden  associated  with  a  complex  geometry 
will  be  organized  internal  to  the  computer  code.  For  example,  the 
operator  need  not  instruct  the  code  that  two  plates  are  attached  to 
form  a  convex  or  concave  structure.  The  code  flags  this  situation 
by  recognizing  that  two  plates  have  a  common  set  of  corners  (i.e., 
a  common  edge).  So  if  the  operator  wishes  to  attach  two  plates  together 
he  needs  only  define  the  two  plates  as  though  they  were  isolated. 
However,  the  two  plates  will  have  two  identical  corners.  All  the  ge¬ 
ometry  information  associated  with  plates  with  common  edges  is  then 
generated  by  the  code.  The  present  code  also  will  allow  a  plate  to 
intersect  another  plate  as  shown  in  Figure  7.  It  is  necessary  that 
the  corners  defining  the  attachment  be  positioned  a  small  amount  through 
the  plate  surface  to  which  it  is  being  connected.  These  details  will 
be  illustrated  in  the  examples  of  Section  VII. 


Figure  7.  Data  format  used  to  define  a  flat  plate  intersecting 

another  flat  plate. 
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In  defining  the  plate  corners  it  is  necessary  to  be  aware  of 
a  subtlety  associated  with  simulating  convex  or  concave  structures  in 
which  two  or  more  plates  are  used  in  the  computation.  This  problem 
results  in  that  each  plate  has  two  sides.  If  the  plates  are  used  to 
simulate  a  closed  or  semi-closed  structure,  then  possibly  only  one 
side  of  the  plate  will  be  illuminated  by  the  antenna.  Consequently, 
the  operator  must  define  the  data  in  such  a  way  that  the  code  can  infer 
which  side  of  the  plate  is  illuminated  by  the  antenna.  This  is  ac¬ 
complished  by  defining  the  plate  according  to  the  right-hand  rule. 

As  one's  fingers  of  the  right  hand  follow  the  edges  of  the  plate  around 
in  the  order  of  their  definition,  his  thumb  should  point  toward  the 
illuminated  region  above  the  plate.  To  illustrate  this  constraint 
associated  with  data  format,  let  us  consider  the  definition  of  a  rec¬ 
tangular  box.  In  this  case,  all  the  plates  of  the  box  must  be  specified 
such  that  they  satisfy  the  right-hand  rule  with  the  thumb  pointing 
outward  as  illustrated  in  Figure  8.  If  this  rule  were  not  satisfied 
for  a  given  plate,  then  the  code  would  assume  that  the  antenna  is  within 
the  box  as  far  as  the  scattering  from  that  plate  is  concerned. 

A  missile  structure  is  being  analyzed  as  shown  in  Figure  6. 

The  far  field  of  this  structure  is  classically  defined  as  a  range  ex- 
ceeding  2(D+2AI)  /X.  If  this  code  is  used  to  compute  the  far  field 
pattern,  the  user  should  supply  the  minimum  far  field  range  requirement. 
This  should  be  done  in  that  excessive  ranges  can  create  numerical  dif¬ 
ficulties  within  the  code.  In  fact,  it  is  suggested  that  the  user 
input  the  range  as  (D+2AI)  /X  for  the  far  field  missile  case.  In  any 
event  if  mysterious  far  field  patterns  result  using  this  code,  one 
should  reduce  the  range  in  order  to  examine  the  numerical  accuracy 
of  the  solution. 

In  the  "PG:"  command,  if  LATACH(MPX)=T  ( i . e . ,  the  plate  is  at¬ 
tached  to  the  fuselage),  the  program  assumes  that  the  first  and  last 
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plate  corners  PVC(N,1,MPX)  and  PVC(N,MCMX,MPX)  are  positioned  on 
the  fuselage.  The  user  must  define  the  geometry  accordingly.  How¬ 
ever,  he  need  not  exactly  attach  the  first  and  last  corner  to  the  fuse¬ 
lage  in  that  the  code  will  extend  the  edges  and  reset  the  first  and 
final  corner  points  on  the  fuselage  as  shown  in  Figure  9. 


Figure  9.  Data  format  used  to  define  a  flat  plate 
intersecting  a  fuselage. 

Using  the  "SG:"  command,  it  is  necessary  that  -90°  <  PHSA(MS)  <_  90°. 
In  case  the  antenna  is  placed  on  the  bottom  part  of  the  fuselage,  the 
user  must  redefine  the  geometry  such  that  4>s( PHSA(MS) )  falls  within 
the  required  angular  range.  As  shown  in  Figure  10,  the  fuselage  can 
be  turned  upside  down  by  changing  AF  into  CF',  BF  into  DF',  etc.  to 
correct  this  situation. 

The  ray  tracing  is  always  a  tedious  problem  in  pattern  analysis. 

It  assumes  that  the  higher-order  diffraction/reflection  fields  from 
the  fuselage  are  small  in  which  case  they  can  be  neglected.  However, 
the  code  automatically  simulates  fuselage  blockage  using  two  additional 
plates  (i.e.,  one  along  the  x-axis  and  the  other  along  the  y-axis). 

Thus,  even  though  higher-order  interactions  between  structures  and 
the  fuselage  are  not  added  in  computation,  their  absence  will  be  apparent 
in  the  results. 
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(a)  (b) 


I  NOTE  :  AF  *  CF1 
BF  =  DF' 

CF  *  AF‘ 

DF-BF'  y 

Figure  10.  The  original  (a)  and  redefined  geometry  (b)  to  adapt 
to  the  data  format  in  defining  the  antenna  location. 

Finally,  it  must  be  kept  in  mind  that  the  antenna  should  be  kept 
at  least  a  quarter-wave length  away  from  any  diffracting  edge.  In  fact 
all  dimensions  should  be  at  least  a  quarter  wavelength. 
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V.  PROGRAM  OUTPUT 

The  basic  output  from  the  computer  code  is  a  line  printer  listing 
of  the  results.  Recall  that  the  results  of  the  program  are  the  E@p 
and  E^p  radiation  pattern  values.  In  order  to  again  describe  these 
pattern  components,  let  us  consider  the  various  principal  plane  pat¬ 
terns  treated  in  the  previous  section.  The  computer  code  allows  for 
a  rotation  of  coordinates  such  that  one  can  take  a  pattern  about  an 
arbitrary  axis.  This  information  is  input  to  the  code  using  the  spherical 
angles  THC,  PHC.  The  geometry  that  applies  for  each  of  the  roll,  ele¬ 
vation,  and  azimuth  patterns  of  the  next  section  is  illustrated  in 
Figure  II.  Note  that  the  0p  and  <j>p  angles  are  defined  relative  to 
the  rotated  pattern  coordinates  and  that  they  change  as  THC,  PHC  is 
changed.  Thus,  Eflp  is  the  theta  component  of  the  field  (i.e.,  E9p=£*9p) 
in  the  pattern  coordinate  system.  Likewise,  E<i>P=?  *<J>p.  The  total  radiated 
electric  field  is  denoted  by  £. 

In  addition  to  the  printed  results,  one  has  the  option  of  ob¬ 
taining  a  set  of  polar  patterns.  If  LPLT=T  in  the  input  data  list, 
the  program  will  automatically  plot  the  E0p  and  E^p  polar  patterns. 

These  patterns  are  plotted  such  that  the  outer  ring  corresponds  to 
the  pattern  maximum  in  each  case.  This  polar  plot  routine  was  used 
to  plot  the  data  presented  in  the  next  section. 
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(b)  ELEVATION  PLANE  COORDINATES  (THC=90°, 

PHC  =  ±  90° ) 


(c)  AZIMUTH  PLANE  COORDINATES  (  THC  =  90°, 

PHC  =  0°  ) 


Figure  11a.  Illustration  of  pattern  coordinates  for  the  principal 
plane  pattern  calculations  in  roll  plane  model. 
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(a)  ELEVATION  PLANE  COORDINATES  (THC  =  0°, 

PHC  s  0°  ) 


(b)  AZIMUTH  PLANE  COORDINATES  (THC=90°, 

PHC  *  0°  ) 


Figure  lib.  Illustration  of  pattern  coordinates  for  the  principal 
plane  pattern  calculations  in  elevation  plane  model 
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VI.  APPLICATION  OF  CODE  TO  SEVERAL  SIMPLE  EXAMPLES 

The  following  examples  are  used  to  illustrate  some  features  of 
the  computer  code.  Various  geometries  and  input  data  are  used  to  ex¬ 
amine  the  effect  of  separate  terms.  Note  that  the  patterns  are  plotted 
in  decibels  with  each  division  being  10  dB,  and  the  labeling  is  not 
included  on  the  plots  in  some  cases. 

Example  1: 

Consider  the  roll  plane  pattern  of  an  axial  slot  antenna  mounted 
on  a  cylinder  as  shown  in  Figure  12. 

X 

z 


Figure  12.  Geometry  of  an  axial  slot  mounted 
on  a  circular  cylinder. 

The  input  data  is  given  by 

UN: 

3 

SG: 

1 

0.  ,0. 

0.08,0.51,  0. ,0.25,1 

l.,0. 
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FG: 

6.75,6.75,6.75,6.75 

0. ,0. ,0. 

PD: 

0. ,0. ,90. 

0,360,1 
600., 12. 

EX: 

The  E^p  roll  plane  pattern  is  plotted  in  Figure  13.  The 
pattern  is  not  plotted  because  it  is  of  negligible  value. 
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Example  2: 

As  shown  in  figure  14,  consider  now  the  case  that  4  intakes  are 
attached  to  the  fuselage  with  an  axial  slot  antenna  mounted  on  it. 

The  input  data  is  given  by 

UN: 

3 

SG: 

1 

1.32,0. 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

0. ,1.32,90. 

0,360,1 

600..  12. 

IG: 

6. 75, 9. 125, 2. 375, 2. 375, 4, T 
EX: 

The  calculated  and  measured  E.D  roll  plane  pattern  are  compared 

<PK 

in  Figure  14.  Note  that  the  source  location  is  adjusted  to  correspond 
to  the  experimental  result. 
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Example  3: 

Following  Example  2,  the  Input  data  for  an  elevation  plane  pattern 
is  given  by 

UN: 

3 

SG: 

1 

1.32,0. 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

90..  91. 32. 90. 

0,360,1 

600..  12. 

IG: 

6. 75, 9. 125, 2. 375, 2. 375,4, T 
EX: 

The  calculated  and  measured  EQp  elevation  patterns  are  compared 
in  Figure  15. 
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Example  4: 

Following  Example  2,  but  with  an  axial  slot  antenna  mounted  on 
the  intake  instead  of  the  fuselage,  as  shown  in  Figure  16,  the  input 
data  is  given  by 

UN: 

3 

SG: 

1 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

0. ,0. ,90* 

0,360,1 

600..  12. 

IG: 

6. 75, 9. 125, 2. 375, 2. 375, 4, F 

EX: 

The  calculated  and  measured  E^p  roll  plane  patterns  are  compared 
in  Figure  16. 
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Example  5: 

Following  Example  4,  the  input  data  for  an  elevation  plane  pat¬ 
tern  is  given  by 

UN: 

3 

SG: 

1 

3.76,0. 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

90..  93. 76. 90. 

0,360,1 

600..  12. 

IG: 

6. 75, 9. 125, 2. 375, 2. 375, 4, F 
EX: 

The  calculated  and  measured  E0p  elevation  plane  patterns  are 
compared  in  Figure  17. 
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Example  6: 

In  order  to  examine  the  effect  of  the  separate  terms  on  the  total 
pattern,  the  "TO:"  command  can  be  used.  Consider  again  the  geometry 
in  Figure  14  but  with  <j>$=0,  the  input  data  for  the  total  and  separate 
field  pattern  is  given  by 

a)  Total  field: 

UN: 

3 

SG: 

1 

0.,0. 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

0. ,0. ,90. 

0,360,1 

600. . 12. 

IG: 

6. 75, 9. 125, 2. 375, 2. 375, 4,T 
EX: 

b)  source  field: 

TO: 

F,F,F 

T,T 

T,F,F,F,F,F,F 

1,6,1 

1,4,1 

1,4,1 

1,4,1 
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1.4,1 

1.4,1 

1.4,1 

UN: 

3 

SG: 

1 

0.  ,0. 

0.08,0.51,0. ,0.25,1 

1.. 0. 

PD: 

0. ,0. ,90. 

0,360,1 

600..  12. 

IG: 

6. 75,9. 125, 2. 375, 2. 375, 4, T 
EX: 


The  fourth  line  can  be  adjusted  to  yield  the  reflected  or  diffracted 
field,  etc.  The  resulting  E  p  roll  plane  pattern  is  plotted  in  Figure 
18. 


(c)  REFLECTED  FIELD  (d)  DIFFRACTED  FIELD 

Figure  18.  Roll  plane  patterns  for  an  axial  slot 
located  at  <t>s=0°. 


NADC -80096 -30 


VII.  APPLICATION  OF  CODE  TO  AIRCRAFT  SIMULATIONS 

A  typical  aircraft  model  contains  a  cylinder  and  flat  plates. 

Two  such  configurations  are  the  elevation  and  roll  plane  models  which 
can  be  used  to  calculate  most  patterns.  The  elevation  plane  model 
consists  of  a  composite  elliptic  cylinder  chosen  to  approximate  the 
actual  aircraft  side-view  drawing.  Aircraft  front-view  drawing  is 
for  roll  plane  model.  The  input  parameters  used  to  describe  these 
two  aircraft  models  are  given  in  Figures  10  and  11.  Although  not  shown 
in  Figures  10  and  11,  one  can  simulate  various  aircraft  structures 
such  as  stores,  pylons,  wing  tanks,  etc.  by  multiple  flat  plates. 

This  code  allows  for  two  different  methods  for  defining  one  plate 
to  be  attached  to  another:  1)  edge  to  edge  attachment  and  2)  edge 
to  surface  attachment.  Edge  to  edge  attachment,  as  illustrated  in 
Figure  19,  often  requires  that  a  plate  edge  be  defined  as  two  or  three 


Figure  19.  Edge  to  edge  plate  attachment. 

colinear  edges  as  the  program  identifies  this  mode  of  attachment  only 
by  finding  two  identical  pairs  of  corners.  Note  that  the  corners  must 
be  consecutive  on  both  plates  which  means  there  actually  exists  an 
edge  between  them.  In  the  case  of  edge  to  surface  attachment,  one 
plate  is  defined  as  penetrating  a  short  distance  through  the  surface 
of  the  second  plate  as  illustrated  in  Figure  20.  The  program  then 
defines  the  new  junction  edge  and  eliminates  the  smaller  portion  of 
plate  #1  behind  plate  #2.  Here  care  must  be  taken  to  assure  that  the 
new  junction  edge  is  completely  contained  within  the  bounds  of  plate 
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Figure  20.  Edge  to  surface  plate  attachment. 


#2,  and  no  where  nearer  than  a  half  wavelength  or  so  to  an  edge  of 
plate  #2.  Now  it  will  be  shown  how  these  two  attachment  methods  are 
'used  to  simulate  aircraft  wings. 

A  roll  plane  model  for  a  Cessna  402B  is  shown  in  Figure  21  with 
the  input  file  given  in  Figure  22.  This  is  an  illustration  of  a  case 
in  which  edge  to  edge  plate  attachment  is  used.  Since  the  engine  extends 
beyond  both  the  forward  and  rear  limits  of  the  wing,  it  is  possible 
to  define  the  four  points  at  which  the  engine  and  wing  intersect. 

Then  for  the  sake  of  simplicity,  the  part  of  the  wing  beneath  the  engine 
is  eliminated,  which  is  permissible  since  the  source  is  mounted  such 
that  it  will  not  directly  illuminate  that  area.  Note  that  the  sides 
of  the  engine,  plates  #2  and  #4  each  have  six  edges. 

A  case  illustrating  edge  to  surface  plate  attachment  is  that 
of  the  Beecher aft  Baron.  The  roll  plane  model  for  this  aircraft  is 
shown  in  Figure  23  and  the  input  file  given  in  Figure  24.  Here  the 
engine  does  not  extend  behind  the  rear  limit  of  the  wing,  necessitating 
the  defining  of  the  rear  part  of  the  engine  to  be  mounted  flush  on 
the  wing.  To  make  such  an  attachment,  the  entire  junction  edge  must 
be  contained  within  the  boundaries  of  the  wing.  Thus  one  must  use 
the  tab  poked  slightly  through  as  shown  in  the  three  dimensional  (3D) 
view.  Note  that  there  is  a  distinct  separation  between  the  front  wing 
edge  and  the  engine  according  to  the  rules  of  edge  to  surface  attachment. 
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UN I  INCHES  CESSNA  4021  LOC  1  ROLL  U/  FUEL  TANKS  PBiLEFT  FUEL  TANK  PLATE  «1 


3 

TOt 

f,f,f 

F,T 

PSl 

0.,0.,?0. 

0,340,1 
4200., 3. 

PBiLEFT  UIN6  (INNER  PART) 
4,T 

-41.8,24.0,0. 

—34.8,48. ,0. 

—34.8,48. , 43 . 

-41. 8, 24. ,43. 

PBiLEFT  ENGINE  (RIGHT  SIDE) 
4,F 

—34 .8,48. ,0. 

-34.8,48. ,-40. 

-21. 8, 48. ,-40. 

-21. 8, 48. ,83. 

—34.8,48. , 85 . 

—34.8,48. ,43. 

PBiLEFT  ENGINE(TOP) 

4»F 

-21. 8, 48. ,-40. 

-21. 8, 104. ,-40. 

-21. 8, 104. ,85. 

-21. 8, 48. ,85. 

PBiLEFT  ENGINE  (LEFT  SIDE) 
4,F 

—34.8,104. ,45. 

—34 .8,1 04 . , 85 . 

-21. 8, 104. ,85. 

-21. 8, 104. ,-40. 

-34.8,104. ,-40. 

-34.8, 104. ,0. 

PBiLEFT  ENGINE  (FRONT) 

4»F 

-21. 8, 48. ,-40. 

-34.8,48. ,-40. 

-34.8, 104. ,-40. 

-21. 8, 104. ,-40. 

PBiLEFT  U1N6(0UTER  PART) 

*,F 

-34.8,104. ,0. 

-24. 8, 213. ,0. 

-24.8,213. ,30. 

-34.8, 104. ,45. 


4,F 

-24.8,213. ,0. 

-14.8,21?. ,-32. 

-14.8,21?. ,82. 

-24.8,213. ,50. 

P6iLEFT  FUEL  TANK  PLATE  82 
4,F 

-14.8,21?. ,-32. 

—8 .8,235. ,0. 

-8.8,235. ,50. 

-14.8,21?. ,82. 

S6t 

1 

0.,-10. 

.414, .828,0. ,.25,3 

1 .  ,0. 

F6: 

8.2,24. ,230. ,24. 

0. ,0. ,0. 

EXt 


Figure  22.  Input  data  for  Cessna  4028. 
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UNl  INCHES 
3 

TO: 

F  ,F,F 
F,T 
PS  l 

0 ■ ,0 « ,80 • 

0,300,1 

4200..  5, 

PBi  LEFT  IIING 
4,T 

-10. ,23.5,-20. 

10. . 227. ,-10. 

10. . 227. .27, 

-10. ,23.5,00. 

P8:  RIGHT  UINS 

4,T 

-10. ,-23.5,00. 

10. , -227. ,27. 

10. , -227. ,-10. 

-10. ,-23.5,-20. 

PSi  LEFT  EN6INE  (RIGHT  SIDE) 
0,F 

-5.0,50.5,-02. 

8.4,50.5,-02. 

8.4.50.5.20. 

-8.0,50.5,20. 

-8. 0,50. 5, -8. 5 
-5. 0,50. 5, -8.5 

P6s  LEFT  ENGINE  (TOP) 

4,F 

8.4.50.5, -02. 

11.4.84.5, -02. 

11.4.84.5.20. 

8.4.50.5.20. 

P6t  LEFT  ENGINE  (LEFT  SIDE) 
0,F 

-5.0,84.5,20. 

11.4.84.5.20. 

11.4,84.5,-02. 

-2.0,84.5,-02. 

-2. 0,84. 5, -7. 

—5.0,84.5,— 7 . 

P6:  LEFT  ENGINE  (FRONT) 

4.F 

-2.0,84.5,-02. 

11.4.84.5, -02. 

8.4.50.5, -02. 

-5.0,50.5,-02. 


P6i  RIGHT  EN81NE  (LEFT  SIIE) 
0,F 

-5. 0,-50. 5, -02. 

—5. 0,-50. 5, -8.5 
—8.0, -50.5, -8.5 
-8.0,-50.5,20. 

8.4, -50.5,20. 

8. 4,  -50. 5, -02.0 

PSi  RIGHT  ENGINE  (TOP) 

4»F 

8. 4,  -50. 5, -02. 

8.4, -50.5,20. 

11.4, -84.5,20. 

11.4, -84.5,-02. 

PSi  RIGHT  ENGINE  (RIGHT  SIDE) 
6,F 

-5.0,-84.5,20. 

—5 .0, -84 .5,-7. 

-2. 0,-84. 5, -7. 

—2 . 0, -84 .5, -02 . 

11.4, -84.5,-02. 

11.4, -84.5,20. 

PSt  RIGHT  ENGINE  (FRONT) 

4,F 

—2 .0, -84 .5 , -02 . 

—5 . 0, -50 .5,-02 . 

8.4,  -50.5,-02. 

11.4, -84.5,-02. 

S6t 

1 

0. ,-77. 

.414, .828,0. ,.25,3 

t.,0. 

F6: 

11.7,17.5,22.0,17.5 

0.,0.,0. 

EX: 


Figure  24.  Input  data  for  Beechcraft  Baron. 
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Example  7: 

Consider  a  crossed-slot  mounted  on  the  fuselage  of  a  KC-135  air¬ 
craft.  The  computer  model  of  the  KC-135  is  illustrated  in  Figures 
25  and  26.  The  input  data  is  given  by: 


UN l  CHANGE  TO  INCHES 
3 

FGiKC-133  AIRCRAFT  STUDY 

3.  ,3. ,3. ,3. 

0.  ,  0 ■  fO. 

SB i  CROSSED-SLOT  FORWARD  OF  UIN6S 
2 

43. . 8.34 

0. ,0.3,0. ,0.25,1 

1.. 0. 

43. . 8.34 

0. ,0.5,90. ,0.23,1 
1 .  ,90. 

P6i  LEFT  UING 
«,T 

-1., 3. ,12.31 
-1., 28, 5, 34. 41 
-1.,  28. 5, 40. 41 
-1., 3. ,24.41 
P6i  RIGHT  UING 

4, T 

— 1 . , "3. ,24.41 
-1., -28. 5,40.41 
—1 • ,-28.5,34.41 
-1. ,-3. ,12. 31 
PD «  ROLL  PATTERN 
0.  ,0, ,90. 

0,340,1 

1000. . 187.3 

EX i  EXECUTE  PROGRAH 
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Figure  25.  Computer  simulated  roll  plane  model  for  a  KC-135.  The 
antenna  is  located  forward  of  wings  (off  center). 
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Figure  26.  Computer  simulated  roll  plane  model  for  a  KC- 135.  The 
antenna  is  located  over  wings  (off  center). 


The  E.„,  E-_  radiation  patterns  in  the  roll  plane  model  for  the  forward 
9P  “P 

and  over  wing  locations  are  given  in  Figure  27  and  Figure  28,  respectively. 
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Figure  27.  Roll  plane  pattern  of  KC 
mounted  in  the  forward  1< 
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Example  8: 

Consider  a  monopole  mounted  on  the  fuselage  of  a  C- 141  aircraft. 
The  computer  model  of  the  C-141  is  illustrated  in  Figure  29.  The  input 
data  is  given  by: 


UN:  CONVERT  TO  FEET 

2 

F6:  C-141  AIRCRAFT  ROLL  PLANE  NODEL 

7.37,7.37,7.37,7.37 

0. ,0. ,0. 

PB:  UIN6 
5,T 

A. 0,7. 37, 11.0 
A. 0,78. A, 45. 44 
A. 0,78.4,55.23 
A. 0,30. 7, 43. 

A. 0,7. 37, 41. 5 
PS:  UING 
5,T 

A. 0,-7.37,41 .5 
A. 0,-30. 7, 43. 

A. 0,-78. A, 55. 2A 
A. 0,-78. A, 45. 44 
A. 0,-7. 37, 11. 

SG:  NONOPOLE 
1 

0..18.4 

0. ,0.5,0. ,0. 25 , 3 

1.,0. 

PD:  AZINUTH  CONICAL  PATTERN 

90.0,0.0,50.0 

0,340,1 

1000.0,2.52 

EX: 
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UNi  CONVERT  TO  FEET 
2 

F6i  C-141  AIRCRAFT  ELEVATION  PLANE  HODEL 
7.37,27.45,7.37,104.5 

o.,o.,o. 

PBi  VERTICAL  STABILIZER 
4,T 

0.-104.45,1.93084 

24.58, -111.21,1. 

24.58, -95.85,0. 

11.47, -84.79,0. 

P6t  VERTICAL  STABILIZER 
4,T 

11. 47,  -84.79/0. 

24.58, -95.85,0. 

24.58, -111.21,-1. 

0., -104. 45, -1.93084 
P6i  T-TAIL 

4,F 

24.58, -111.21,1. 

24.58, -114.42,25.3 

24.58, -110.45,25.3 

24.58, -95.85,0. 

P6i  T-TAIL 

4,F 

24.58, -95.85,0. 

24.58, -110.45,-25.3 

24.58, -114.42,-25.3 

24.58, -111.21,-1. 

PDt  ELEVATION  ANB  THETA  CLOSE  TO  90  CONICAL  PATTERN 

90. . 0. .70. 

0,340,1 

1000. . 2.52 
S6l  HONOPOLE 
1 

-44. ,0. 

0.2, 0.5,0. ,0.25,3 

1.. 0. 

DO:  BOUBLE  TERNS 
4 

3,3,3, 1 
3, 3, 3, 4 

4. 1.4. 3 

4. 1 .4.4 

EX l  EXECUTE  PR06RAN 


The  elevation  plane  model  is  used  to  get  azimuthal  conical  patterns 
(6p=60°,  70°,  80°,  90°)  and  elevation  pattern  (T-tail  effect). 

Note  that  the  roll  plane  model  is  used  to  compute  the  azimuthal  conical  patterns 
(6p=10°,  20°,  30°,  40°,  50°,  100°,  110°).  The  radiation  patterns  are 
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Figure  29.  Computer  simulated  model  for  the  C-141  aircraft, 
(a)  roll  plane  model,  (b)  elevation  plane  model, 

(c)  roll  plane  model. 


compared  with  measured  results  made  on  a  -|q  scale  model  of  a  C-141 
at  General  Dynamics  (San  Diego,  California)  in  Figures  30,  31,  32  and 
33.  The  T-tail  effect  on  the  elevation  pattern  is  illustrated  in 
Figure  34. 
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Figure  30.  Scale  model  (1/10)  of  C-141  aircraft  at 
General  Dynamics  (San  Diego,  California) 
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Figure  31.  Comparison  of  measured  (dash  line)  and  calculated 
(solid  line)  results  of  azimuthal  conical  pattern  at 
(a)  yl0°,  (b)  ^=20°,  (c)  y30d,  (d)  ^=40°. 
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Figure  32.  Comparison  of  measured  (dash  line)  and  calculated 
(solid  line!  results  of  azimuthal  conical  pattern. 

(a)  9p=50cr,  (b)  ep=60°,  (c)  0p=7O  ,  (d)  0p=8O°. 
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Figure  33.  Comparison  of  measured  (dash  line)  and  calculated 
(solid  line)  results  of  azimuthal  conical  Dattern. 

(a)  ep=90°,  (b)  ep=ioo  ,  (c)  0p=l 10  , 

(d)  elevation  pattern. 
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Figure  34.  Calculated  elevation  pattern  (a)  without  T-tail  and 
vertical  stabilizer  (dash  line),  without  T-tail  but  with 
vertical  stabilizer  (solid  line),  (b)  with  T-tail 
and  vertical  stabilizer. 
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Example  9: 


Consider  a  monopole  mounted  on  the  bottom  fuselage  of  a  F-4 

aircraft.  The  computer  model  of  the  F-4  is  illustrated  in  Figure  35. 

Note  that  the  pattern  coordinate  system  (with  subscript  p)  has  been 

turned  upside  down  on  purpose.  If  the  antenna  is  mounted  on  the  belly 

of  an  aircraft,  then  one  needs  to  turn  the  aircraft  upside  down.  This 

is  done  to  insure  that  one  defines  the  geometry  associated  with  the 

bottom  of  the  aircraft  different  from  the  top,  in  that  most  aircraft, 

especially  military  aircraft,  have  a  definite  shape  change  from  top 

to  bottom.  In  that  the  F-4  was  turned  upside  down,  there  will  exist 

a  coordinate  system  difference  between  the  true  aircraft  and  our  model. 

In  fact  for  the  ep=75°  pattern  computed  here,  the  actual  pattern  angle 

from  the  vertical  is  180°-75°  or  105°.  With  this  difference  in  mind 

the  input  data  is  given  by: 

UN i  INCHES  F-4  AIRCRAFT 
3 

TO* 

F,F,F 

F,T 

PD:  AZINUTH  CONICAL  PATTERN 
90.0,0.0,73. 

0,340,1 

3000.. . 373 

F6:  F-4  AIRCRAFT  ROLL  PLANE  MODEL 

20. . 20. .20. .20. 

0.,0.,0. 

SB:  N0N0P0LE 
1 

0.,0. 

.414, .828,0. ,.25,3 

t.,0. 

P8i  ATTACHED  PLATE 
4,T 

19., -?. ,-150. 

30. , -13. ,-150. 

40., -4.,-1 50 • 

40 « , 4 . , —150. 

30. . 13. ,-130. 

19. . 9. ,-150. 

PB:  ATTACHED  PLATE 
4»T 
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13. . 18. ,-344. 

13..  230., -375. 

13. . 230. ,-327. 

13. . 50. ,-138. 

13. . 50. ,-75. 

13..  18,j— 75. 

P8l  ATTACHED  PLATE 
4,T 

13. , -18. ,-75. 

13.,  —50. , — 75 . 

13. , -50. ,-138. 

13. , -230. ,-327. 

13. , -230. ,-375. 

13. , -18. ,-344. 

P8<  ATTACHED  PLATE 
4,T 

13. . 18. ,-75. 

13. . 50. ,-75. 

-37. ,50. ,-75. 

-37. ,18. ,-75. 

P6:  ATTACHED  PLATE 
4»T 

-37. ,-18. ,-75. 

—37 ■ ,—50. ,—75. 

13. , — 50. ,—75. 

13. , -18. ,-75. 

P6t 

4,F 

41. . 72., -PI. 

15. . 72., -91. 

15. . 72. ,-143. 

41. . 72. ,-143. 

Pfli 

4,F 

41. , -72. ,-143. 

15. , -72. ,-143. 

15. , -72. ,-91. 

41. , -72. ,-91. 

P6i 

4,F 

51. . 127. ,-131. 

15. . 127. ,-131. 

15. . 127. ,-354. 

91.  ,127. ,-354. 

P6i 

4,F 

31. , -127. ,-354. 

15. , -127. ,-354. 

15. , -127. ,-131. 

51. , -127. ,-131. 

EXi 
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Figure  35.  Computer  simulated  model  for  the  F-4  aircraft 

(roll  plane  only). 


Figures  36,  37,  38,  39  and  40  indicate  the  E0p  pattern  contributions 
associated  with  different  portions  of  the  aircraft  structure.  One 
may  gain  a  better  feeling  for  the  model  through  using  this  kind  of 
process  in  complex  shape  pattern  computations.  The  final  computed 
E@p  pattern  is  compared  with  a  measured  result  obtained  at  the  RADC 
Newport  antenna  range  in  Figure  41. 

The  F-4  is  obviously  a  complex  airframe  with  many  structures 
attached  to  the  basic  frame.  If  an  antenna  system  is  to  be  developed 
for  such  an  aircraft,  the  antenna  designer  needs  to  know  how  the  various 
structures  associated  with  the  airframe  affect  the  resulting  patterns. 
Simply  using  measured  data  this  is  not  always  available  in  that  a  detailed 
complete  model  is  usually  experimentally  examined.  However  using  this 
code,  the  antenna  designer  can  examine  the  effects  of  various  structures 
as  done  here  for  the  F-4  and  later  for  the  A-10. 
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Figure  36.  Calculated  F-4  azimuthal  conical  pattern  for 
*■  cylinder  only  model  at  0p=lO5°  frequency  =  375  MHz. 
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Figure  38.  Calculated  F-4  azimuthal  conical  pattern  for 
cylinder,  two  wings  and  blockage  underneath  the 
fuselage  model  at  e  =105  ,  frequency  =  375  MHz. 
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Example  10: 


Consider  four  monopoles  mounted  on  the  belly  of  an  A- 10  aircraft. 
The  computer  model  of  the  A-10  is  illustrated  in  Figure  42.  The  input 
data  is  given  by: 


TOi 

F,F,F 

FrT 

UN:  INCHES  A-10  AIRCRAFT 
3 

POi  AZIMUTH  CONICAL  PATTERN 

90.0,0.0,75.0 

0,360,1 

6000.0.17.576 

P6i  ATTACHED  PLATE 

4,T 

0.0,0.425,-4.62 

0.0,1.58,-4.62 

0.0,1.58,-2.69 

0.0,0.425,-2.69 

PBi  ATTACHED  PLATE 

M 

0. ,-0.425,-2.69 
0., -1.58, -2. 69 
0., -1.58, -4. 62 
0. ,-0.425,-4.62 
P6i 
4,F 

0.54,1.58,-4.62 

0.54,1.58,-2.69 

0.0,1.58,-2.69 

0.0,1.58,-4.62 

PBi 

4,F 

0.54,1.58,-4.62 

0.54,1.86,-4.62 

0.54,1.86,-2.69 

0.54,1.58,-2.69 

PBi 

4,F 

0.54,1.86,-2.69 

0.54,1.86,-4.62 

0.0,1.86,-4.62 

0.0,1.86,-2.69 
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Ml 

4,F 

0.0,1.86,-4.42 

-0.3,5. 0,-4. 34 

-0.3, 5. 4, -3. 11 

0.0,1.80,-2.04 

Mi 

3,F 

0.18,1.713,-2.02 

0.0,1.58,-2.04 

0.34.1.38,-2.04 

Ml 

3,F 

0.18,1.715,-2.02 

0.34,1.38,-2.04 

0.34,1.80,-2.04 

Ml 

3,F 

0.18,1.715,-2.02 

0.0,1.80,-2.04 

0.0,1.38,-2.04 

P8t 

3, F 

0.18,1.713,-2.02 

0.34,1.80,-2.04 

0.0,1.80,-2.04 

P8t 

4, F 

0.0,1.38,-2.04 

0.0,1.80,-2.04 

0.0,1.80,-4.02 

0.0,1.38,-4.02 

P6i 

4,F 

0.0,-1. 38,-4.02 

0.0,-1. 38,-2.04 

0.54,-1.38,-2.04 

0.54,-1.58,-4.02 

PBi 

4,F 

0.54,-1.38,-2.04 

0.34,-1.80,-2.04 

0.34,-1.80,-4.02 

0.34,-1.38,-4.02 

PBi 

4,F 

0.0, -1.80,-2.04 
0.0, -1.80,-4.02 
0.34, -1.80, -4.02 
0.34,-1.80,-2.04 
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Hi 

4,F 

0.0, -1.84, -2.49 
-0.5, -5. 4, -3. 11 
-0.3, -3.4, -4.39 
0.0,-!. 84,-4.42 
P6i 
3,F 

0.18,-1.713,-2.02 

0.34,-1.38,-2.49 

0.0, -1.58, -2.49 

P6i 

3,F 

0.18,-1.713,-2.02 

0.34,-1.84,-2.49 

0.54,-1.38,-2.49 

P6i 

3,F 

0.18,-1.713,-2.02 
0.0, -1.58, -2.49 
0.0, -1.84, -2.49 
P8i 

3. F 

0.18,-1.713,-2.02 
0.0, -1.84,-2.49 
0.34,-1.84,-2.49 
P8i 

4, F 

0.0, -1.38,-4.42 
0.0, -1.84,-4.42 
0.0,-1 .84,-2.49 
0.0,-1 .38,-2.49 

38 t  H0N0P0LE  ARRAY  WITH  COUPLING  EFFECT  INCLUBE1  ON  RELATIVE  CURRENT  VALUE 
4 

0.0,-0.824 

0.0,0.0,0.0,0.23,3 

0.272,14.0 

0.0, -1.14 

0.0,0.0,0.0,0.25,3 

1.0, -3.0 

73.47,-1.328 

0., 0.0, 0.0, 0.25, 3 

0.272,14.0 

-73.47,-1.328 

0.0,0.0,0.0,0.25,3 

0.272,14.0 

F8l  A-10  AIRCRAFT  ROLL  PLANE  NOBEL 

0.117,0.425,3.07,0.425 

0.0, 0.0, 0.0 

EXi 


[ 
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Figure  42.  Computer  simulated  model  for  the  A-10  aircraft 
(roll  plane  model  only). 


Note  that  four  monopoles  are  spaced  a  half  wavelength  apart.  Mutual 
coupling  between  the  antennas  Is  significant  for  this  spacing  and  cannot 
be  neglected  In  the  pattern  computations.  This  problem  can  be  solved 
for  this  geometry  using  a  thin  wire  moment  method  code  to  treat  the 
four  closely  spaced  loaded-dipoles  (50  ohms  at  center  of  each  dipole). 

The  relative  current  distributions  on  the  four  monopoles  are  approximately 
the  same  as  dipoles  In  free  space  If  curvature  of  the  surface  where 
monopoles  are  mounted  Is  smooth  (Image  theory  then  holds  here  for  engine¬ 
ering  purposes).  Figures  43,  44  and  45  Indicate  the  E0p  pattern  contri¬ 
bution  from  different  portions  of  the  aircraft  structure.  The  final 

computed  Ea_  pattern  Is  compared  with  a  measured  result  obtained  at 
up 

the  RADC  Newport  site  In  Figure  46. 
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Figure  43.  Calculated  A-10  azimuthal  conical  pattern  for  cylinder 
only  model  at  ep*105  ,  frequency  ■  300  MHz. 
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Figure  44.  Calculated  A-10  azimuthal  conical  pattern  for  cylinder 
and  two  flat  plates  model  at  6  *105°, 
frequency  *  300  MHz.  p 
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Figure  46.  Comparison  of  measured  (solid  line)  and  calculated 

/  _ i  s  j _ \  _ .1  _ i _ i.i _ i _ i _ 1 _ _ 


i)  final  results  of  azimuthal  con 
at  9p=105  ,  frequency  *  300  MHz. 
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